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Nonlinear, electrocatalytic swimming in the presence of salt 

Benedikt Sabass-^ and Udo Seifert-^ 

//. Institut fiir Theoretische Physik, Universitdt Stuttgart, 70550 Stuttgart, Germany 

A small, bimetallic particle in a hydrogen peroxide solution can propel itself by means of an electrocatalytic 
reaction. The swimming is driven by a flux of ions around the particle. We model this process for the presence 
of a monovalent salt, where reaction-driven proton currents induce salt ion currents. A theory for thin diffuse 
layers is employed, which yields nonlinear, coupled transport equations. The boundary conditions include a 
compact Stern layer of adsorbed ions. Electrochemical processes on the particle surface are modeled with a 
first order reaction of the Butler- Volmer type. The equations are solved numerically for the swimming speed. 
An analytical approximation is derived under the assumption that the decomposition of hydrogen peroxide 
occurs mainly without inducing an electric current. We find that the swimming speed increases linearly with 
hydrogen peroxide concentration for small concentrations. The influence of ion diffusion on the reaction rate 
can lead to a concave shape of the function of speed vs. hydrogen peroxide concentration. The compact 
layer of ions on the particle diminishes the reaction rate and consequently reduces the speed. Our results are 
consistent with published experimental data. 

PACS numbers: 47.63.mf, 47.61.-k, 82.45.-h, 87.16.Uv 
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I. INTRODUCTION 

Motion of micrometer-sized swimmers, directly driven 
by an inhomogeneous chemical surface reaction, has re- 
cently attracted much scientific interest. While the mech- 
anism was originally suggested as a mode of biological 
self propulsion^^^, currently studied swimmers are mostly 
artificial. A variety of systems featuring this chemical 
self-propulsion, have been investigated^'^. The employed 
reactions often involve hydrogen peroxide which decom- 
poses on the metal^^ or enzyme-coated^^ surfaces of the 
swimmers. Other mechanisms are, e.g., oxidation and 
reduction of glucoseii, hydrolyzationi^ or bromination^^ 
at droplet surfaces, and thermally induced phase separa- 
tion around a swimmer—. Finally, a number of theoret- 
ical studies have been conducted with generic chemical 
reaction scheme^—. In previous theoretical worl^ii^ 
we have investigated the efficiency of this surface-driven 
propulsion since its exciting potential for use with micro- 
machines has been amply demonstrated^ —. Chemically 
driven microswimmers can also be used as model systems 
to investigate the physical principles of non-equilibrium 
diffusive motion^^^^. Mutual interaction^— and the 
influence of confinement on the swimming^"— are a fur- 
ther active field of research. 

A number of physical driving mechanisms for hydrogen 
peroxide based propulsion have been suggested. One of 
them is diffusiophoresis, which is particle motion through 
its interaction with a concentration gradien t ""^^i^^i^^ . An 
alternative, that is seen to dominate for large-scale 
swimmers, is the recoil of oxygen bubbles^i^. Pres- 
sure waves were also suggested as a possible origin of 
motior^. Finally, it has been found that metallic 
microswimmers swim through self-electrophoresis^^^^^. 
Self-electrophoresis is the swimming of a charged par- 
ticle through electric interaction with a self-generated, 
charged environment. 

The work in this article is motivated by the ex- 
periments with metallic swimmers done by Sen, Mal- 



louk, Wang and co- worker s^SiiU following their work with 
Paxton et al^^. There, bimetallic microrods (often a 
platinum-gold particle) were immersed in an hydrogen 
peroxide {H2O2) solution. The surface of the rods cat- 
alyzes a decomposition of H2O2 into oxygen O2 {g) and 
water. It has been found that the mechanism for elec- 
trokinetic decomposition of H2O2 involves an electric 
current inside the microrod, which causes a concentra- 
tion gradient of ions in the fluid surrounding the swim- 
mer. Since the rod interacts through an electric potential 
with the ions, a gradient of the latter leads to swim- 
ming. The swimming speed increases linearly with the 
concentration of H2O2 until it saturates above concentra- 
tions of about 5 %. An interpretation of this saturation 
has been given in the framework of a Michaelis-Menten- 
like surface kinetics. Wang et al.^^ conducted a compre- 
hensive study using different combinations of metals for 
bimetallic microswimmers. A clear correlation between 
the mixed potential difference between the two metals 
and the swimming speed was found. By choice of opti- 
mal materials^ for the anode, cathode and intermediate 
part of the swimmer, the speed could even be increased to 
more than 50 /im/s. In order to understand the electroki- 
netic driving mechanism quantitatively, a few theoretical 
studies have been conducted. Sundararajan et al.— com- 
pared experimental data with a numerical model where 
the rate of the cation exchange was fixed. Yariv— con- 
ducted a thorough study of the case of a thin diffuse 
layer of charges surrounding the particle. His linear re- 
sponse theory included surface reaction kinetics explic- 
itly. Moran and Posner^^ investigated a full numerical 
model for rather thick, diffuse layers of ions around the 
swimmer. Most recently, they included a surface reaction 
that is second order in concentration^^. This work 
provided a closed, quantitative model that could be di- 
rectly compared to experimental data. Their results in- 
dicated a quadratic increase of the swimming speed with 
the concentration of H202-, which, however, is not ob- 
served in the experiment. Building on the work described 
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above, we here investigate a complementary, non-linear 
model for a charged double layer structure that is much 
thinner than the radius of the swimmer. 

An estimate of the ion concentrations^^ around the 
swimmers suggests that these systems are operating quite 
far away from equilibrium. On the other hand, the lin- 
ear speed vs. concentration relation at low and moderate 
H2O2 concentrations does hint that a theory containing 
linear response elements may still be valid. Therefore, 
we choose a combination of numerical modeling and ana- 
lytical approximation to capture the non-equilibrium be- 
havior of the system and provide simple explanations for 
qualitative trends. 

Any model of self-electrophoretic swimming is ham- 
pered by a general lack of knowledge concerning the com- 
plex details of the surface reaction mechanism. In par- 
ticular, the dependence of surface charge accumulation 
on the reaction rate is largely unknown. Since the elec- 
tric potential of the swimmer directly influences its swim- 
ming speed, a better understanding of the effect of charge 
screening and ion adsorption is highly desirable. To pave 
the way for further investigations of these problems we 
model the system as simple as possible, demonstrating 
only the generic aspects. 

II. THE MODEL 

A. Swimmer in an ionic solution 

We employ the commonly used hydrodynamic model 
for electrophoretic effects^^. Throughout this publica- 
tion, we will designate dimensional variables with a tilde 
( ~ ). Dimensional constants are written in calligraphic 
letters. The swimmer, a spherical particle with radius 7^, 
is placed in an infinitely large container (Fig. [1]). Swim- 
ming speed will be denoted by The fluid surrounding 
the particle is assumed to be incompressible and Newto- 
nian with a constant viscosity 77. Mass flow velocity is 
denoted by v. The fluid contains ions that carry only 
one unit charge each. The concentrations of salt anions 
and cations in the solution are denoted by c*~ and c*+, 
respectively. The hydroxide {0H~) and proton {H^) 
concentrations, resulting from the spontaneous dissocia- 
tion of water and from the electrocatalytic process, are 
denoted by and c^+, respectively. All types of ions 
are assumed to have the same diffusion constant V. Con- 
centrations of cations and anions are denoted summarily 
by 



The ions couple to an electric field <j). In the bulk, far 
away from the swimmer, we have charge neutrality and 
the concentrations of ions are equal c~(oo) = c+(oo). 
The swimmer is axially symmetric. Therefore we use a 
spherical coordinate system aligned in the direction. 



where r is the distance from the particle center and is 
the inclination angle. The unit vectors of the spherical 
system are denoted by e^, e^. 

We non-dimensionalize length with the swimmer ra- 
dius IZ and the concentrations of ions with the bulk 
concentration of salt ions C = c*'+(oo) = P'~(oo). 
Energies are normalized with the thermal energy scale 
kbT. Accordingly, the electric potential (j) is made non- 
dimensional with the thermal voltage k^T/Ze 25 mV. 
Throughout the publication we set Z = 1. Diffusive 
fluxes are non-dimensionalized with VC/IZ. The velocity 
scale of the flow v and the swimming speed U is given 
by a typical speed U. The Peclet number associated with 
transport of cations and anions is given by Pe = UIZ/V. 




FIG. 1. Schematic representation of the model for a self- 
electrophoretic swimmer with radius IZ moving with swim- 
ming speed U . The electric potential outside the swimmer 
is denoted by (j). The inset shows an enlarged version of our 
double layer model. The swimmer surface is covered by an 
immobile, compact layer of salt ions, termed Stern layer. The 
inner Helmholtz plane (iHp) and outer Helmholtz plane (oHp) 
delimit the modeled region of the Stern layer with thickness 
As. Bordering to the outer Helmholtz plane we have the dif- 
fuse layer, where dilute solutes form a charge cloud around 
the swimmer. The diffuse layer has a thickness A. 



B. Poisson-Boltzmann and diffusion equations 

The fluxes of cations and anions are given in dimen- 
sionless form by 

■^{H,i},+ = _^c^H,^},+ _ + + pg^^{ff,i},+ ^ (2) 

Diffusion equations guarantee the conservation of sub- 
stance in the fluid. They determine the concentration 
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field and read 



V.jW},+ = 0, 



(4) 
(5) 



The Poisson-Boltzmann equation provides a connection 
between tiie ion density and ttie electric field. It reads in 
dimensionless form 



(6) 



The lengthscale of the potential (j) in Eq. © is given by 
the dimensionless Debye length A, which we define as 



A 



ekbT 1 



(7) 



with A/a being the Avogadro constant and e the permit- 
tivity of the solution. 



where the parameter Xg is the dimensionless thickness of 
the Stern layer between inner and outer Helmholtz plane. 

The electrocatalytic process leading to a decomposi- 
tion of H2O2 is supposed to take place at the inner 
Helmholtz plane. It leads to an absorption and emis- 
sion of ions with a rate non-dimensionalized 
with the flux CV/IZ. Boundary conditions, determining 
the concentrations of ions, result as 



e,r'+(l,i9) = a(^), 
e,r±(l,i?) = e,j^'-(l,i?)=0 
(r ^ 00,^) = 1, 



c^'±(r^ 00,1?) = (5. 



(11) 
(12) 
(13) 
(14) 



The constant 5 is the relative concentration of ions 
in the bulk. It is defined by 



(oo)/C. 



(15) 



C. Stem layer and boundary conditions 

In the bulk, far away from the swimmer, the electric 
potential takes on an uniform value which we set to zero 



(j){r ^ 00,1?) = 0. 



(8) 



The presence of salt in the solution suggests some form 
of adsorption of ions onto the metal surface. The inset 
of Fig. [1] illustrates the employed model of a Stern (com- 
pact) -layer of immobile ions covering the metal surface^^. 
All the charge in the Stern layer is assumed to be concen- 
trated in one plane, called inner Helmholtz plane (iHp). 
This plane almost coincides with the metal surface. We 
do not consider the details of the structure beyond the 
inner Helmholtz plane. The boundary between the Stern 
layer and the diffuse layer is termed outer Helmholtz 
plane (oHp). The voltage difference between inner and 
outer Helmholtz plane is called Stern layer voltage drop 
A05(i9). The electric potential at the inner Helmholtz 
plane V equals A^5(i?) plus the change of the potential 
between outer Helmholtz plane and infinity 0(1, '^9): 



y = A,^,(t?) + (/.(f,t?). 



(9) 



Since we have an electric current flowing through the 
swimmer, the potential V may not be strictly constant. 
However, the electric conductivity of the metal swimmer 
is much higher than the conductivity of the surrounding 
solution. We therefore assume that the variation of V is 
negligible. A priori, the potential change over the Stern 
layer, A(j)s{'d), is unknown. If it depends on the potential 
right outside the immobile layer, it becomes i^-dependent. 
Since one postulates that the Stern layer outside the in- 
ner Helmholtz plane is charge free. Gauss's law can be 
employed to arrive at 



A(ps{if) = -As — \r=l, 



(10) 



In order to keep the charge of the swimmer constant the 
exchange of cations at its surface must satisfy 



e,j^'+(l,i?) dAr^i 



0. 



(16) 



D. Hydrodynamic equations 

The mass flow velocity v is determined by the Stokes 
equation for overdamped motion 

V • cr = V^v -\/p = g{c-^ - C-) (17) 

where p is the pressure and a = Vv + (Vv)^ — pi is the 
hydrodynamic stress tensor. We have also defined the 
non-dimensional constant 



nchT 

3 = T^' 

The boundary conditions on the fluid flow are 

v(r ^ oo,^^) = U, 
v(0,^) = 0. 



(18) 



(19) 
(20) 



where the swimming velocity U is calculated from the 
balance of forces on the particle 







= J ^/(c+ -C-) V0dV + J aerdAr=i. (21) 



For a body force g (c+ — c~) Vcj) that is independent of 
the fluid velocity v one can employ Teubner's formula!^ 
to calculate the swimming speed as 



U ■■ 



9_ 



(22) 
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E. Redox reactions at the surface of the swimmer 

The electrochemical reactions leading to the decom- 
position of H2O2 at the swimmer's surface are compli- 
cated, especially since a number of inhibitory processes 
are involved. Experimental evidence supports the view 
that the reaction rate depends non-linear ly on the H2O2 
concentrations^, most visibly for CH2O2 ^ ^ The in- 
ferred (Michaelis-Menten-like) kinetics for a microswim- 
mer agrees with other findings^^. This kinetics is proba- 
bly due to saturation of active sites at the surface. Here 
we do not model the saturation but focus on low concen- 
trations of educts and products. Experiments with com- 
paratively large (mm-sized) electrodes suggest that the 
decomposition of hydrogen-peroxide at metal surfaces in 
acidic medium happens via a formation of oxides at the 
surfaces^. As also discussed in App. [Aj this intermedi- 
ate may serve to justify an effective description in terms 
of a first order surface reaction. Fig. [2] illustrates that 
protons combine with H2O2 for the reduction of H2O2' 
Therefore, the rate of consumption is assumed to 
be ~ c^'~^{l,'d)cH202' Electrocatalytic oxidation is sup- 
posed to supply the protons with a rate proportional to 
CH2O2' 



oxidation reduction 




FIG. 2. Sketch of the electrokinetic process taking place at 
the swimmer's surface as described in Sec. Ill El Note that the 
swimmer is found to move with the oxidizing end forwards. 



satisfy detailed balance. We choose rates ~ i , 
where the transfer coefficient, modifying A^^ in the ex- 
ponents, is set to 1/2. For simplicity, we neglect further 
direct effects of the compact salt layer, that could, e.g., 
reduce the catalytic activity by covering the active sites. 
Evidently, (j) also modifies the reaction rate by inffuenc- 
ing the concentration of protons right outside the Stern 
layer c^'+(l, "d). This dependence is however implicit and 
it does not appear in the rate equation. The overall rate 
of proton exchange at the swimmer's surface takes the 
form 

a(^) = - (i^i + ki cos^) e-^^c^'+(l, i?)ci^202 

+ (i^2 - ^2 COS ^) e 2 CH2O21 

with Ki > ki and K2 > k2. The first term in Eq. (|23|) 
models the catalytic reduction of H2O2 while the second 
term models the provision of ions through an oxida- 
tion. The constants ki and /c2 refiect the difference of 
redox potentials between the two metals composing the 
swimmer. Note that the linear dependence of a{'d) on 
the H2O2 concentration in the bulk merely rescales all 
rate constants in the same way. The bulk concentration 
of H2O2 is a model parameter. We choose to express 
CH2O2 (dimensionless) fractions of 1% wt/v in order 
to facilitate the comparison with experimental data. In 
not modeling the H2O2 concentration explicitly we have 
assumed that a diffusive supply of this reagent is not the 
bottleneck of the reaction mechanism and that the en- 
counter of H2O2 with the surface can be incorporated in 
the effective rate constants. 



III. THEORY FOR THIN DIFFUSE LAYERS 

If the Debye length is small, one can attempt an ex- 
pansion of the concentration and potential fields in or- 
ders of A, which largely simplifies the nonlinear equa- 
tions. We follow here the classical approach outlined, 
e.g, by Hunter— which was used recently in the group 
of Bazan t^^i^^ to model time dependent electrochemical 
processes. 

A. The nonlinear problem 



The redox reactions at the swimmer's surface cause 
proton fiuxes and involve a charge transfer. There- 
fore, the electric potential infiuences the reaction rate. 
We model this dependence with the established Butler- 
Volmer kinetics^—. It is assumed that the reaction 
takes place at the inner Helmholtz plane, almost coin- 
ciding with the metal surface of the swimmer (see Fig. 
[1]). Since the Stern layer is only described in an effective 
way, the energetic cost of passing from the outer to the 
inner Helmholtz plane must be incorporated explicitly 
into the reaction rates. This is commonly done by intro- 
ducing rates that depend exponentially on A(j)s{i^) and 



The technical details of the nonlinear double layer 
model are lengthy. We therefore defer them to App. 
[B1 Here we only sketch the mathematical procedure and 
present the result for the swimming speed. The field 
variables are expanded in powers of A as 

c''± =c^'± + Aci'± + ..., 
c^.±=Co^'±+Acf± + ..., (24) 

^ = (/)0 + A (/)1 + . . . . 

In the following, we only consider the terms of lowest or- 
der, O(A^). The strategy is now to solve the non- linear 
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field equations both in an outer region, far away from the 
surface, and inside the diffuse layer. The inner and outer 
solutions are then matched to each other, such that an 
uniformly valid solution can be constructed. We distin- 
guish the outer variables from the inner ones through a 
hat ( on the former. Also, the potential (f)o is written 
as sum of an inner and outer solution 



(25) 



The quasi one-dimensional equations for the inner vari- 
ables can be solved analytically while the outer equations 
must be solved numerically. The result for the electric 
potential in the diffuse layer is 
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(26) 



where we defined y = (r — 1)/A. The function ^{1^) is 
determined by the boundary conditions on the swimmer 
surface (see App. [Bj). B{i^) is proportional to the overall 
concentration of ions at the outer boundary of the diffuse 
layer 

B{^) = {cl^'^M + c!^'-{l,^) + c};+{l,^) + (l,^))/2. 

(27) 

Given (numerical) solutions for 7/^0 and the potential at 
the swimmer surface Vq, we can use Eq. (|9]) to calculate 
the potential drop across the diffuse layer A(j)o{i^) as 



reversed with the sign of the zeta potential. Therefore, 
the first term can be interpreted as the electrophoretic 
part. The second term does not change sign with the 
potential A(/)o since -41n (cosh (A0o/4)) = 21n(l -7^) 
is negative semi-definite. Its dependence on a gradient 
in overall concentration suggests an interpretation of the 
second term as diffusiophoretic contribution. 



B. Analytical approximation 

We assume now that the redox potential difference 
between the two sides of the swimmer is considerably 
smaller than its overall redox potential. This case can be 
modeled by setting Ki^ K2 ki^ k2. In a steady state, 
the constant H2O2 decomposition will modify the elec- 
tric charge of the swimmer. If, additionally, ki and /c2 
are neglected, no net cation flux can take place since a 
radially symmetric flux violates the charge conservation. 
The concentration of protons near the surface then fol- 
lows a Boltzmann distribution ~ ^exp [— 0o,o(^)], where 
the second index 0, here at ^0,0, indicates that we are 
considering the solution of order (/c^^, /C2). In this limit, 
Eq. ([23|) becomes 



= 



-Kie ^(5e-^°'°(°) 



CH2O2 



K2e 2 CH2O21 



(28) 



where the voltage drop across the Stern layer A(j)s{'d) is 
calculated through Eq. (p!Q]) . A0o('^) can be interpreted 
as an angle-dependent zeta potential. The final result for 
the swimming speed reads, in terms of the above defined 
variables, 



(30) 

It follows that the swimmer's potential in steady state 
Vo,o = A0,,o + ^0,0(0) = In {K15/K2) (31) 

does not depend explicitly on the H2O2 concentration to 
lowest order. The swimming speed is calculated in App. 
O We find 



U^gX 



2 2 I n (1 + 70) 
1 + S 



I 



ao,i cosi^dcosi^ 



(32) 



JO 



d-d 



41n(cosh(^))^i^sin^(.)]d.. 

(29) 

We will omit the index at Uq in the following since we 
do not deal with higher orders in A. Remarkably, Eq. 
(f29|) has the same form as the celebrated formula by 
Dukhin and Derjagui n^^i^^ , which was derived for the 
linear response regime. The only difference is that the 
potential drop across the diffuse layer can not be ap- 
proximated by the equilibrium zeta potential and is now 
7^- dependent. The dependence of U on A^o allows to sep- 
arate electrophoretic and diffusiophoretic contributions 
in the lowest order swimming speed —. The first term in 
the integrand of Eq. ([29|) changes sign with A^o? similar 
to electrophoresis, where the particle swimming is also 



where the integral is a measure of the dipole moment 
of the first order emission/absorption rate ao,i- We 
also defined 70 = tanh(^o,o(0)/4). Eq. (j32j) is similar 
to previously derived linear response formulas—. When 
inserting ao,i, which is calculated in App. [Cl we have 



;ln(l + 7o) ^fc + feiff) CH2 



1 



28 i 



K2CH2O2 



(33) 



For a weak potential outside the compact layer, (/)o,o ^ 1, 
one can estimate the dependence of the swimming speed 
on the particle potential through ln(l + 7o) — {Vq — 
A(/>.,o)/4. 



IV. NUMERICAL VALUES FOR THE MODEL 
CONSTANTS 

The thermal enersv scale is hT = 4.11 x 10"2^J. For 
the solution viscosity, we choose 77 = 9 x 10~^Pas. For 
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the permittivity of the solution around the swimmer, we 
choose 6 = SOcyak- Typicahy, the experimental swim- 
mers have lengths in the 300 nm - 5 /im range. We set 
the radius of our swimmer to 

7^ = 1 /im. (34) 

A velocity scale can be defined, which is independent of 
the concentration scale. We set 

UgX^ = p^c,20.Q^m/s. (35) 
87r?7 

In this work we rely on the addition of extra salt to the 
solution to explicitly justify the usage of the thin diffuse 
layer model. Where not explicitely mentioned otherwise, 
we set the bulk concentration of salt ions in the following 
to 

C = 5 X 10"^ mol/L. (36) 

For a monovalent salt, e.g., NaNOs as used by Paxton 
et al.— , according to Eq. ([7|), we have a thickness of 
the diffuse layer of 12 nm. Since this is much smaller 
than 1Z = 1 /im our theory for A <C 1 is expected to work 
well. The absolute thickness of the Stern layer is usually 
assumed to be in the order of a molecular diameter—, 
which might be about 0.3 nm. 

For propulsion of our microswimmers we consider a 
[lO~^ ... 5] % solution of H2O2 in water. Percent of 
H2O2 concentration are given in wt/v (1 % wt/v cor- 
responds to lOg H2O2 per liter solvent, molar mass of 
H2O2'' 34.02 g/mol). The pH value decreases here al- 
most linearly with an increase of H2O2 concentration^. 
It lies roughly between 7 and 5. This sets the concen- 
tration of bulk protons and hydroxide ions to c^'^(oo) 
[lO~^ . . . 10~^] mol/L. The relative concentration of 
ions in the bulk is then in the range 8 = c^(oo)/C 

[io-^..io-i]. 

For the diffusion constant, we choose V = 7 x 
10~^m/s. The ion fluxes at the surface are then non- 
dimensionalized with 

CV/n = 3.5 X 10""^ mol/m^ s. (37) 

A measured^^ evolution rate of O2 for a 3.7 % solution of 
H2O2 is 8.7 X 10~^mol/m^ s. Within our reaction model, 
Eq. ([23|) . this rate roughly determines the magnitude of 
the second term, responsible for the oxidation of H2O2' 
For CH2O2 — 1 (1 %) we estimate that K2 CH2O2 is on the 
order of 10~^. We assume that its value scales linearly 
in the range of employed H2O2 solutions. The constants 
Ki^ /ci , k2 are chosen such that the we obtain reasonable 
values for the swimming speed and the proton fluxes. 

V. RESULTS 

In order to facilitate a comparison with experimen- 
tal data, we often display dimensional quantities in the 



following. The dimensional speed is calculated as = 
UU = 20.6 /im/s X U/(gA^). Results for the analytical 
approximation of U (Eq. ([33])) are also shown where ap- 
plicable. 

A. Reaction induced concentration distortion 

Fig. [3] shows exemplary concentration fields of posi- 
tive ions around the swimmer. While the concentrations 
of in Figs. [3] a,c are directly determined through 
the surface reactions, the inhomogeneous distribution of 
positive salt ions seen in Figs. [3] b,d is an indirect ef- 
fect. The positive salt ions accumulate oppositely to 
the ions. They present a non- negligible contribu- 
tion to the local charge balance. In Figs. [3] a,b we have 
(/ci, /C2) <C (i^i, ^^2), i-e, the decomposition of H2O2 
happens fairly homogeneously around the swimmer. In 
this case the concentration profile is antisymmetric and 
follows the distribution of surface reactivity. In Figs. [3] 
c,d we have set ki = Ki and /c2 = i^2, i-e., the reaction 
happens mainly in an asymmetric way, where oxidation 
and reduction take place on different sides of the swim- 
mer. Also, the concentration of ions in the bulk is 
not assumed to be excessive. Then the concentration 
variation around the swimmer feeds back on the overall 
reaction rate, Eq. ([23|) . Consequently, the antisymmet- 
ric distribution of the cations seen in Figs. [3]a,b is lost. 
This reaction induced concentration distortion is a direct 
consequence of the usage of the mass action law to model 
the surface reaction. In our previous work we have found 
that it is important for the calculation of the swimming 
speedi^. 

B. Swimming and H2O2 concentration, fixed pH 

Fig. m demonstrates the dependence of the swimming 
speed U on the concentration of H2O2 with a fixed pH 
value. A fixed pH value means that the availability of 
in the bulk, determined by ^, does not change with i^202 
concentration. Eq. ([33|) approximates the numerically 
calculated swimming speed well for ki < Ki and /c2 < 
K2. However, the condition k2 < Ki for the validity of 
the analytical approximation seems here unnecessary. 

The curves in Fig. (jH) have a concave shape, which 
is not attributed to a saturation of the catalytic sur- 
face in our model. Rather, it results from a limita- 
tion of the reaction rate through the diffusion around 
the swimmer. Quantitative understanding of this ef- 
fect can be gained from the approximative Eq. (|33|) . 
Here we find a leveling off of the swimming speed since 
the rate of ion fiux through the swimmer saturates for 

K2CH2O2 ^ 2(5 e 2 . This estimate is seen to hold also 
for the results from the nonlinear numerics. A negative 
Stern layer voltage drop A(j)s shifts the saturation of the 
speed to higher values of ch^o^i while also reducing the 
speed (see Sec. IV Ep . Figs. [5l [6] show how the average 
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FIG. 3. Colorcoded concentrations of ions and positive 
salt ions. Figs. a,b and c,d each have have the same color 
coding. CH^o^ = 1 (1%), Ki = 0.0185, K2 = 0.185, 5 = 0.1, 
and As = 0. For Figs. a,b: ki =0.1/^1, /c2 = 0.1 K2 resulting 
in Vq = -4.61 and U = -0.0398. For Figs. c,d: /ci = Ki, 
k2 = K2 resulting in Vq = -5.03 and U = -0.485. 



FIG. 4. Swimming speed U vs. H2O2 concentration at fixed 
pH value. Ki = 0.005, 3 — 0.1, As = 0. Symbols are nu- 
merical results. Full lines are approximative speeds from Eq. 
(|33l). (Red box): K2 = 0.1, /cs = K2, ki = Ki. (Open 
red box): K2 = O.l, k2 = 0.25^2, ki = 0.25 i^i. (Green 
diamond): K2 = 0.05, k2 = K2, ki = Ki. (Open green 
diamond): K2 = 0.05, /c2 = 0.25 i^^2, ki - 0.25 ii^i. (Blue tri- 
angle): K2 = 0.025, k2 = K2, ki = Ki. (Open blue triangle): 
K2 = 0.025, k2 = 0.25/^2, ki = 0.25 ii^i. 



reaction rate {\a{'d)\) and the electric potential on the 
swimmer's surface Vq change with the increase of H2O2' 
Note that (|a(i9)|) can not be identified with the measur- 
able oxygen evolution rate. {\a{'d)\) vanishes for equal 
redox potentials on both sides of the swimmer, i.e., if 
ki = k2 =0, while the oxygen evolution rate in this limit 
is, within our model, an unknown constant. 



C. Swimming and H2O2 concentration, variable pH 

Since the concentration of H2O2 also influences the 
pH value we plot in Fig. [7] the speed vs. CH2O2 with 
i^202 -dependent S. The pH value decreases almost lin- 
early with an increase of H2O2 concentratior^. There- 
fore we set 6 = Sq -\- Si CH2O2 where the dissociated 
water molecules at a pH value of 7 contribute Sq = 
(10~'^mol/L)/C. The slope of the function is roughly 
estimated as Si = 0.02. The reaction rate is not ex- 
pected to saturate here as seen from Eqns. (|33IC20p . The 
speed becomes an almost linear function of the H2O2 con- 
centration. A remaining concave tendency of the func- 
tions is due to a reduction of the Debye length and a 
change of the particle potential with the pH value. The 
curves in Fig. [71 are very similar to what is measured 
experimentally^!^. 




HjOg [%] 

FIG. 5. Average ion current (|Q;('i^)|) vs. H2O2 concentration 
at fixed pH value. Symbols are the same as in Fig. [H 



D. Swimming and H2O2 reduction rate 

In Fig. Owe investigate the role of the rate constants 
of H2O2 reduction, Ki and /ci, for the swimming speed. 
We keep the quotient of Ki and ki constant and Eq. 
([33]) therefore suggests that the swimming speed is only 
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FIG. 6. Electric potential of the swimmer Vb vs. H2O2 con- 
centration at fixed pH value. Symbols are the same as in Fig. 
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FIG. 8. Dependence of swimming speed U on the rate con- 
stant of H2O2 reduction Ki at CH2O2 = 1 (1 %) with k2 = K2, 
ki = Ki, and A, = 0. (Black box): S ^ 0.1, K2 ^ 0.05. 
(Blue Dot): 6 = 0.1, K2 ^ 0.025. (Violet circle): S = 0.01, 
K2 = 0.025. (Red triangle): 6 = 0.1, K2 = 0.01. 
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FIG. 7. Dependence of swimming speed U on H2O2 con- 
centration where the pH value varies with ch202- We use 
6 = 0.002 + 0.02 CH202- Symbols are the same as in Fig. [H 



influenced through its dependence on the particle poten- 
tial Vb via the function 7. The approximation for the 
swimmer's potential, Eq. (|3T]) . in turn suggests that Vq 
depends logarithmically on K1S/K2. The numerical data 
qualitatively supports a scaling of U with this logarithm. 
The magnitude of the swimming speed is seen to decrease 
with an increase of Ki. The negative sign of U is in ac- 
cordance with the experimental finding that the swim- 
ming occurs in direction of the i^202 -oxidizing end^. 
However, for KiS > K2 we also find a reversion of the 
swimming direction, which results from the sign change 



of K). 



E. Effect of the Stern layer 

The Stern layer modifies the swimming speed for Vb < 
by reducing the reaction rate and the potential change 
in the diffuse layer, appearing in Eq. (j29j). Fig. [9] demon- 
strates that U decreases for our rate constants almost lin- 
early with the average Stern layer voltage drop {A(j)s{i^)). 
The approximation, Eq. (|33|) . does not yield satisfactory 
absolute values for U in Fig. [9]since we have here ki = Ki 
and /c2 = K2. However, Eq. (|33|) can be expanded for 
small A(j)s. It then predicts the slopes of the decrease of 
U with A(j)s fairly well. The inset of Fig. [9] demonstrates 
that A(j)s depends non-linearly on on the thickness of the 
layer A^/A. This non-linear relationship occurs since the 
Stern layer voltage drop is determined through a tran- 
scendental equation (see appendix, Eq. ()B30p ). 

The swimming speed scales within the thin diffuse 
layer model, Eq. (|29|) , as the inverse of the overall bulk 
ion concentration 1/(C -\- CS). Experimental data sup- 
ports this rough scaling, but also indicates deviations 
from it^^. One reason for a deviation may be the Stern 
layer, which influences the reaction rate. In our simple 
model, the connection between the Stern layer and the 
reaction is given through the appearance of A^^ in the 
Butler- Volmer Eq. (|23]) . In a tempting first approxi- 
mation, we here keep the absolute thickness of the Stern 
layer constant. The variations of the bulk salt concentra- 
tion still affects the Stern layer voltage drop A^^ since 
Ag/A ~ VC. Fig. (fTQ|) demonstrates that this suffices 
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(A(^^) [kj/e] 



FIG. 9. Dependence of swimming speed U on the average 
Stern layer voltage drop (A(/)s('i^)). Inset: dependence of 
(A(/)s('i^)) on relative Stern layer thickness As/ A. CH2O2 — 1 
(1%) with k2 = K2, ki = Ki. (Black box): 6 = 0.1, 
Ki = 0.001, K2 = 0.1. (Blue Dot): S = 0.1, Ki = 0.001, 
= 0.05. (Green star): 6 = 0.1, Ki = 0.001, i^^2 = 0.01. 
(Red triangle): 6 = 0.01, Ki = 10"^; K2 = 0.01. (Dotted 
lines): slopes predicted by an expansion of Eq. (|33|) for small 




FIG. 10. Deviation of U from the standard thin diffuse layer 
scaling with the inverse ion concentration ^ 1/{C -\-C6). Bulk 
salt concentrations are varied as C = n x (5 x 10~^ mol/L) 
with n = 1 ... 10. Speed is multiplied with n{l-\-S) to remove 
the standard scaling. CH2O2 = 3.7 (3.7%), Ki = ki = 0.005, 
K2 = k2 = 0.025/n, 6 = 0.1/n. (Triangle): Xs/X = 0. (Star): 
As/A = 0.005 yn. (Dot): Xs/X = 0.025 v^. (Diamond): 
Xs/X = 0.1 y/a. 



to cause a reduction of the swimming speed compared 
with the standard scaling. A precise measurement of 
this effect, along with a more refined model of the ion 
adsorption is called for to improve the quantitative un- 
derstanding here. 

VI. DISCUSSION 

Effects of salt on self-electrophoretic swimming become 
important when (bio-) technological applications are en- 
visaged since physiologically or otherwise relevant en- 
vironments are often ionic. One of these effects is the 
couphng of salt ion fluxes to the reaction-driven proton 
fluxes. In general, salt ions can not be described by a con- 
stant concentration background since this approximation 
causes non-negligible errors in the swimming speed. A 
full mathematical description leads to a system of non- 
linear equations, that must be solved numerically. We 
chose to focus on the thin diffuse layer approximation 
rather than solving the whole system of equations. The 
advantages of this approach are that the hyodrodynam- 
ics become relatively simple and that a quite demand- 
ingly high density of grid nodes near the surface can be 
avoided. The thin diffuse layer model may describe the 
swimming well, even in absence of electrolytes. Sponta- 
neous dissociation of H2O2 in the bulk can be an impor- 
tant factor here. A reasonable pH value of 6 translates, 
e.g., to a double layer thickness, Eq. ([7|), of about 90 nm. 



This number is an order of magnitude smaller than IZ. 
It justifies the approximative use of the thin diffuse layer 
model for H2O2 concentrations, say, above the 2 % range. 
A contamination of the water with atmospheric carbon 
dioxide has also been found to be important for elec- 
trophoresis in a salt-free environment^^. The somewhat 
unknown details of the reaction mechanism at the surface 
of the swimmer could further contribute to making the 
thin diffuse layer model valid. It is, e.g., possible that 
the reactions at the swimmer's surface lead, beyond our 
model, to a further accumulation of other types of ions 
and radicals. 

A new feature of our model is that we use first order 
Butler- Volmer kinetics for the electrocatalytic decompo- 
sition of H2O2' This reaction leads to a linear increase of 
swimming speed at low H2O2 concentrations. The trans- 
port of protons around the swimmer can limit the reac- 
tion rate and therefore cause a saturation of the swim- 
ming speed at high H2O2 concentrations. We found that 
an expected decrease of the pH value with an increase 
of the H2O2 concentration can widely suppress the lev- 
eling off of the speed since the bulk ions take part in 
the catalytic reduction. Any other mechanism providing 
an excessive amount of ions would equally reduce 
the effect of proton diffusion and lead to a linear speed- 
[H2O2] relation. The influence of proton transport on the 
reaction is not in conflict with an additional limitation 
of the rates though saturation of the catalytic surface in 
experimental systems. 

The Butler- Volmer equation includes a further effect. 
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that is also related to the presence of salt: the Stern layer 
on the surface of the swimmer. In this article, we have 
undertaken a very first step towards an understanding 
of its influence on the swimming speed. The Stern layer 
lowers in our model the swimming speed by modifying 
the potential drop across the diffuse layer and by reducing 
the reaction rate. It is also responsible for a decrease of 
speed with increased salt concentration that is stronger 
than the inverse of the bulk ion concentration. 

The swimmer's potential V plays an important role for 
the swimming speed since it is directly related to the po- 
tential drop across the diffuse layer (the zeta potential). 
For negligible Stern layers, V could possibly be deter- 
mined experimentally by measuring the additional drift 
of actively swimming particles in an external field. In our 
model, we found V to be in most cases much lower than 
— 40mV, which is about the zeta potential of a metal 
particle without the presence of H2O2' This value was 
used previously as a rough estimate^. Such a high poten- 
tial V can only be found within our model if the overall 
oxidation and reduction rates are similar (see Eq. (|3T]) ). 

In this article, we have made a few simplifying assump- 
tions. First, we consider a spherical swimmer with a reac- 
tivity varying like a cosine across the surface. This geom- 
etry is chosen for its numerical robustness since no sharp 
edges and sudden changes are present. On the other 
hand, the experimentally studied swimmers are short, 
bimetallic rods where the reactivity varies quite suddenly 
across the boundary from one metal to the other. A 
second simplificating assumption is that all solutes have 
similar diffusivities. We guess that these simplifications 
together introduce a deviation from the experimental ref- 
erence data by a factor of 1 — 5. 

We have not modeled the adsorption of salt into the 
Stern layer explicitly but fixed its width at a molecu- 
lar lengthscale. The influence ion adsorption on elec- 
trophoretic mobility is a long standing problem. A num- 
ber of different models have been suggested to quantita- 
tively account for different kinds of salt ions^^'^^. The 
common assumption that the Stern layer is indepen- 
dent of salt concentration is somewhat justified for non- 
reactive, metallic interfaces'^ but remains an assumption 
for our system. Its validity may be particularly ques- 
tioned once a variation of the salt concentration changes 
the reaction rate and thereby the swimmer's potential. 
Including an explicit adsorption mechanism is challeng- 
ing. It would be desirable to have more experimental 
data concerning the interplay of salt and H2O2 decompo- 
sition at microswimmers before attempting such a model. 

It is well known that double layers can exhibit lateral 
ion transport if they are highly charged^. If contribu- 
tions of the order \e~^^^^'^ are not negligible, correc- 
tions to the theory must incorporate electro-osmosis 
and electromigration in the diffuse layer. We have ne- 
glected the corrections due to these surface conduction 
phenomena in order to focus on the important trends of 
our model. An investigation of the role of surface trans- 
port for self-electrophoresis would be an important exten- 



sion. A somewhat related topic is the possible occurrence 
of hydrodynamic slip on the surface of the swimmer. This 
slip might serve to enhance the swimming speed^i^. Fi- 
nally, the assumption of dilute solutes could be relaxed. 
However, non-electric interactions and the competition 
for catalytic binding sites would make an analytical de- 
scription of solutes fairly complicated. Molecular dynam- 
ics simulations could offer a promising alternative here. 

Taken together, our model quantitatively supports an 
explanation of the particle swimming in terms of self- 
electrophoresis^SdS and points to a number of interesting 
phenomena governing the details of this intriguing mode 
of micro-motion. 
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Appendix A: Details of the redox reactions 

In this appendix, we illustrate the reaction rate postu- 
lated in Sec. Ill El with concrete, but simplified prototyp- 
ical reaction pathways. We emphasize that the schemes 
discussed below merely present one possible way to jus- 
tify the mathematical model we discuss in the main part. 
Due to the general lack of knowledge concerning the sur- 
face reactions, we can not claim that that these pathways 
are ultimately the dominant ones. However, preliminary 
studies with other rate equations made us confident that 
many qualitative trends presented above are quite robust. 

During the reduction of i^202, ions in the vicin- 
ity of the metal surface combine with electrons form the 
metal substrate S. A possible scheme for this process is 

S + H2O2 + + e- ^ S{OH) + H2O, 

S{OH) + + e- ^ 5 + H2O, 

which is has been used previously for a simulation of 
the reduction on platinum electrode^. Additional pro- 
cesses, that possibly slow down the reduction, are ne- 
glected here. We assume that the concentration of H2O2 
is low, such that the availability of catalyst is not lim- 
iting. Applying the steady state condition, we find that 
the concentration of S{OH) is independent of [H^]. This 
implies that the effective rate equation for consump- 
tion is first order in [H^] and [H2O2]' A second order 
reaction rate can not be excluded, but we deem the sug- 
gested presence of a multi-step decomposition sufficient 
to allow the simpler modeling of the reduction as a first 
order process. 

The oxidation of H2O2 at metal electrodes has been 
studied in great detail by Hall et a L^^i^^ . It is also found 
to depend on surface oxide films. A possible scheme for 
the reduction of the surface film is 

S{0H)2 + H2O2 ^ S{0H)2 ■ H2O2 ^ 5 + 2H2O + O2. 

(A2) 
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Hydrogen ions are released when the surface is oxidized region 



again 



S + 2H2O S{OH)2 + 2iJ+ + 2e- 



(A3) 



Neglecting side processes, one can again assume that the 
reaction is first order in [H202]- Four electron processes, 
as suggested by Wang et al.^^ are, for simplicity, not 
considered explicitly. 



Appendix B: Nonlinear theory for A ^ 

This appendix details the non-linear mathematical the- 
ory of a thin diffuse layer which is used in the main part 
of the article. We define new variables for the concentra- 
tions and expand them in powers of A 



= tuq ^ -\- Ami 

= M(j^'^^+AMp 
= 00 + A 01 + . . . . 



{H,i} 



(Bl) 



Outer variables carry a hat variables for the fields 
inside the diffuse layer do not carry a hat. The potential 
00 is written as sum of an inner and outer solution 



00 = V^o + V^o- 



(B2) 



The equations determining the outer variables Mq 



{H,i} 



and '^o are written in terms of r and 1^. The 



Poisson-Boltzmann equation ([6]) in the outer region reads 



with m 



m 



H 



(^^0 + AV^i • • •) = -j^ (^0 + Ami + A^m2 . . .) , 

(B3) 

which is singular for A ^ and therefore we require 
mo = mi =0. This provides the condition of charge 
neutrality outside the diffuse layer. The number of lowest 
order outer variables is hence reduced by one 



m. 



H 



(B4) 



Note that the condition of charge neutrality does not 
require that the charges of the different ions balance in- 



dividually. The relation m 



H _ 



only holds in the 



bulk, far away from the particle. Concerning the singu- 
lar behavior of Eq. ()B3p we also mention that a charge 
of (vanishing) magnitude ~ A^ ~ 1/C results in a non- 
vanishing lowest order electric field in the outer region. 
Therefore, the electric potential can not be neglected in 
the diffusion equations, Eqns. (|4l5p . outside the diffuse 
layer. Eqns. (|4l5p read to lowest order for the outer 



V • (Vm^ + Mo^VV^o) = 0, (B5) 
V . (- Vm^ + M(5V^o) = 0, (B6) 

V • (vMo^ + V^o) = 0, (B7) 
V • (VM(5 - m^VV^o) = 0, (B8) 

The convection terms do not appear in lowest order be- 
cause the fluid velocity is O(A^) (see below). The bound- 
ary conditions for r ^ 00 are given by the bulk concen- 
tration of ions as 



M^{r ^ 00,1^) = 2 J, 
M(5(r ^ 00,^) = 2, 
m^(r ^ 00,1!}) = 0, 
l^po{r ^ (X),^) = 0. 



(B9) 
(BIO) 
(Bll) 
(B12) 



Near the surface we employ a stretched, radial coordinate 
y = (r — 1)/A. The diffusion equations in the inner region 
are expanded up to O(A^) and read 

= -dy [dyml''^'^ + M^^'^'^dy^o) , (B13) 
= -dy [dyM^""^'^ + mi^'^'^dyi^o) . (B14) 

The lowest order ion flux is 0(1). However, the fluxes 
pertaining to concentrations of O(A^) vanish, since the 
radial derivative introduces a factor of 1/ A 

- (dymj,"-'^ + M^"-'^dy^o) = 0, (B15) 



0. 



(B16) 



The leading order result for the variables ttIq^'^^ {yj'O), 
Mq^'^^ {y,{)) in the diffuse layer are then given by 

(BIT) 
(B18) 

with yet unknown constants. Matching be- 

tween the diffuse layer and the outer region as 
limy^^{mi"''\M^"^'^} = lim.^l{mp^M^^•'>} 
yields 

a^^'^>(i?) -6^^''>(i?) = m^^''^(l,i?), (B19) 
J^'^>(^)+6^^'^>(^) =M(j^'^^(l,i9). (B20) 

The number of constants a^^'^J", in our solutions 

is larger than the number of equations. One needs to 
consider the matching of the 0(A) solutions in order to 
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obtain the full lowest order solution. The result of this 
further matching process are conservation equations for 
the radial solute fluxes, which we do not derive here for 
brevity 

V / r=l 

(B21) 

(^-drM},{r,^)+m^{r,^)drMr,^))^_^=0, (B22) 

(B23) 

(drml^ir, - M^(r, ^)drMr, ^)) = 0. (B24) 

\ / r—1 

In the inner region, the Poisson-Boltzmann equation is 
non-singular. It reads up to O(A^) 



mo 



B{^)smhMy.^). (B25) 



with 



B 



(^)= (M(f(l,i?) + M^(l,i?))/2. 



(B26) 



This equation, determining the non-equilibrium concen- 
tration near the surface, has the same form as the cor- 
responding equilibrium equation. The result for the po- 
tential in the diffuse layer is 



2 In 



.1 -7(t?)e-V^f_ 



(B27) 



The value of the potential at the outer periphery of the 
diffuse layer ipo{l,'d) is determined through the solution 
of Eqns. (|B5IIB8|) . Given i/;o(l, t^), the overall change in 
potential over the diffuse layer A0o(i?) is according to 
Eq. © 

AM'&) = MO,^) = Vo- AM'&) - M^,^)- (B28) 

The voltage drop over the Stern layer is calculated 
through Eq. (p!Q|) . which is in terms of inner variables 



\s dipo{y,^)i 



A 



A 



dy 



\y=o 



(B29) 



Insertion of Eq. (|B29|) into Eq. ()B28|) yields a transcen- 
dental equation for the function 7(1?) 



A 



2 In 



1 + 7(^?) 



1 - 7(t?) 



(B30) 

which now allows the full determination of the potential 
in the diffuse layer. In order to calculate the lowest order 
swimming speed in e^-direction from the above results 



we expand Eq. 
is 



for small A. The lowest order result 



Uo=9^'f I '-^^^[y'dyMy,^)cos^ (B31) 

—ydr&(l)Q ^) sin d] sin ^ d'ddy. 



On inserting Eqns. ()B17|B18|26p including the matched 
constants and using partial integration we flnd Eq. ([29|) 
for the swimming speed. An explicit calculation of U re- 
quires the knowledge of the outer fields at r = 1. There- 
fore, the pertaining differential Eqns. (|B5IIB8p need to 
be solved with the boundary conditions Eqns. (|B9IIB12( 
IB21IIB24]) . Also, the potential of the swimmer, Vq^ de- 
pends through the charge conservation Eq. ([16]) on a 
and must be determined simultaneously with the concen- 
tration fields. The whole task can only be done numeri- 
cally. To this end, we implement an efficient pseudospec- 
tral method^^ where the discretized system of differential 
equations is solved with the Newton- Raphson method. In 
radial direction, we use a Chebychev grid with 30 nodes. 
The grid in angular direction, for '^^ = . . . tt, is uniformly 
spaced with 40 nodes. 



Appendix C: Analytical approximation for U 

In this appendix we provide the details for the ap- 
proximate swimming speed, which is given in Sec. IIIIBl 
Due to radial symmetry and the boundary conditions 
at infinity no ion fiuxes arise in the lowest order where 
ki ~ 0, /c2 ^0. Although being a non-equilibrium steady 
state, the concentrations of ions follow an equilibrium 
distribution throughout the whole system. Solutions for 
the outer field variables, far away from the swimmer, are 
easily determined from Eqns. (|B5IIB8[) with ?/^o,o = to 
be Mq^o 2^, M(5 2, m^Q mj, q 0. The indices 
denote the and k2 order. The inner fields of lowest 
order obey the differential equations 



( 



0,0 ^y9o,Q I 1 



(CI) 
(C2) 



and the solutions are given by Eqns. (|B17IB18p with 



^0,0 



^0,0 



(5, and Q = Q = 1. The potential 



^0,0 



IS 



k,o{y) = 21n 



7oe" 



7oe" 



/B^y 



(C3) 



with Bo = S -\- 1 and 70 = tanh((/)o,o(0)/4). 

For the outer variables of first order in /ci, k2 we de- 



mand niQi = 
read 



-m. 



H 

0,17 



as in App. [Bl The field equations 



V^M(fi = V^M(5^i = 0, 
V'm^i = V'^0,1 = 0. 



(C4) 
(C5) 
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The corresponding boundary conditions are 



'^0,1 ' 



tH 





= <^0,1 


-drMl^\r=l 


= 0, 


- 2(5ar^/'0,l) |r=l 


= Q^o,i 


- 29^^0,1) '■=1 


= 


— 00 — '^Q^\\r—oo 


= 0, 


ll^\r— 00 


= 0. 



(C6) 



Taken together, these equations yield 



Mil = 0, 



Mo^i = ^A„P„(cos^?)r-"-\ 



n=0 



i'o,i=M"/{2S + 2), 



(C7) 
(C8) 

(C9) 
(CIO) 



where the the coefficients of the expansion in Legendre 
polynomials Pn{cosi^) are calculated from 

-S,M(fi|,^i = ^(n + l)A,Pn(cos^) = ao,i. (Cll) 

n=0 

For the corresponding inner variables of first order in 
ki , /c2 we have the differential equations 

= - (S.m^f + M^fdy^o,o + M^^i'dy^o,i) , (C12) 
= - (dyM^f + mQ '^dy(l)o,o + ^o,ody^o,i^ , (C13) 

which result in 

H,i _ H 
"0,1 — ^0,: 

rHA HA-(hr,r, , iH. 



<:l = <fe-*°-« - 6^fe*°-« - 4>o,iM"i\ (C14) 



Mo,i = <1 + e*«'° - 0o,i"^m ■ (C15) 

Matching now as limy^ao M^f = M(^f|r=i, 



we find 6g 1 = — ag 1 = 0, fen'i = and 



«M = X!^"^"(^°^^)■ 



n=0 



This yields for the concentration at the outer edge 
of the compact layer 



^0,1 



(a^i- J(/)o,i)e-^°'°. (C17) 



Having obtained the concentration of cations at the sur- 
face of the swimmer, we are now in the position to cal- 
culate the reaction rate ao,i in first order of /ci, /c2. We 
insert ^ (5e"^0'0 + c^^^ in Eq. (|23|) and employ Eq. 
(|3Q|) for the lowest order reaction rate. The boundary 



condition for the potential Eq. (jSTJ) in order k^, k2 is 
also employed along with the corresponding also in first 
order Stern layer boundary condition 



(CI8) 



to yield 



ao,i e 2 



+ K2Vo,i - cos^ lk2 + k^j^] 

(C19) 



This expression along with Eqns. (|C1H IC16[) allows to 
determine a^^i. Since the swimmer can not emit a net 
flux of cations, we demand that the radially symmetric 
part of ceo, 1 vanishes. Therefore we find with Eqns. (|C1H 
ICTgI) . Ao = 0, 1/0,1 =0, and 



-^0,1 



- fe2 + fclt 

2(5 e ^ +i^2 



(C20) 



In order to calculate the swimming speed we first consider 
the body force in the diffuse layer 



moV^o ~mo,oV(/)o,o + (^0,1 + ^0,1) 9r(/>0,0 
+ (l + (5)(e-^°'° -e^°'°)V^o,i. 



(C21) 



The first, lowest order term is radially symmetric and 
therefore irrelevant for the swimming. On inserting the 
expressions for mo,i, we find 



moV^ ~mo,oV0o,o + a^ie "^^'^ 9r(/)o,o 

+ (l + (5)V((e-^°'°-e^°'°)(/)o,i). 



(C22) 



Since the last term in the body force is a gradient it only 
modifies the hydrodynamic pressure and does cause mo- 
tion of the swimmer in the incompressible limit. There- 
fore, the first order potential (/)o,i is irrelevant for the 
swimming speed. Inserting Eq. (|C22[) with Eq. (|C2Qp 
into Eq. (|B31[) . we find the formula for the swimming 
speed up to first order in the asymmetry of surface reac- 
tivity given in the main text. 
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